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ABSTRA CT 

The theoretical solution for the flow beneath V/STOL aircraft was 
extended to include tilted jet configurations. A laboratory model was 
constructed to test the effect of variation of the parameters governing 
the flow. Free streamline plots, pressure coefficients on the ground 
and fuselage and velocity profiles in the nozzles were determined from 
hot-wire anemometer traverses and micromanometer readings. Experimental 
data compared favorably with the theoretical determinations. 
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a constant 

pressure coefficient 
flow region in the £ plane 

complete elliptic integral of the first kind 

defined by K* (k) = K(k f ) 

modulus of elliptic integral 
/ 2 

defined as \ 1 - k 
pressure 

ambient pressure outside jets 

an angle used to evaluate Q functions 

Reynolds number, AV^/v 

complex variable in T-plane, £ + irj 

complex variable in t-plane 

velocity in x-direction 

ve locity 

velocity at nozzle exit 
average velocity at nozzle exit 
velocity along free streamlines 
velocity in y-direction 
complex potential, § + i\[f 
coordinates in physical plane 
complex variable, x + iy 
tilt angle of nozzle 
plate angle 

flow region in the T-plane 
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width of nozzle 

asymptotic widths of streams flowing to left and to right 

complex conjugate velocity, u-iv 

Theta-function 

function defined in equation (4) 
coordinates of T-plane 
density 

ve locity potential 
stream function 

function defined in Appendix A 



Subscripts : 

A,G,H refer to stagnation points 

r refers to boundary point 
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I. INTRODU TION 



In recent years crowded airport facilities, the search for rapid 
urban transportation and the need for military aircraft capable of 
operating from advanced bases has resulted in a great deal of attention 
being given to V/STOL (vertical/short take-off and landing) aircraft. 

Much research has been performed on various models and full-scale versions 
of these aircraft. Some operational types have been developed. 

When operating near the ground these aircraft are quite similar in 
principle to ground effect machines. Ground effect can be either beneficial 
or detrimental depending upon the configuration of the aircraft. Land 
erosion, damage to objects near the operating area and recirculation of 
debris and hot exhaust gases beneath the fuselage are some of the prob- 
lems involved with hovering near the ground. The reingestion of exhaust 
gases reduces the efficiency of the engines. Recirculating debris may 
also enter the engines or severely damage the underside of the aircraft. 

Some pilots have experienced handling difficulties, drastic loss of power 
and visibility problems as a result of ground effect. 

The recirculating flow underneath the aircraft, however, can also 
produce a lifting force which can significantly reduce the thrust necessary 
to operate the aircraft. Small modifications to the aircraft such as 
tilting of engines can often reduce the unfavorable effects to such an 
extent that the additional lift provided by the recirculating flow can 
produce improved performance. A mathematical model of the flow field 
beneath V/STOL aircraft would permit designers to estimate the effects 
of various modifications analytically. 
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Unfavorable ground effect 





FIGURE 1 
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II. OBJECT 



The theoretical solution for the flow field beneath V/STOL aircraft 
with two vertical two-dimensional jets has been presented* in a series of 
three NASA reports. The purpose of this research was to obtain experi- 
mental verification of the present theory and to extend the theory to 
include other possible configurations. 
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III. THEORY 



An analytical solution to the flow beneath a V/STOL aircraft was 
obtr ined by Goldstein and Siegel [Ref 3]. They considered a two- 
dimensional model of an aircraft with a slot jet exhaust or fan pod in 
each wing. The analysis was performed assuming an incompressible, 
inviscid, isothermal, irrotat iona 1 , steady flow. A flat plate at an 
angle to the ground represented the fuselage. The flow field was mapped 
into a hodograph plane and a complex potential plane. A brief outline of 
this analysis is given in Appendix A. 

So that entrainment could be ignored this solution assumed that the 
nozzles were within a few nozzle widths of each other. For the same 
reason the nozzles were also required to be within a few nozzle widths 
of the ground. The nozzles were assumed to be perpendicular to the 
ground. Since tilting of nozzles is one possible solution to recircu- 
lation problems, a more general analysis should allow the nozzles to be 
at some angle to the perpendicular. 

Considering the flow field shown in Figure 2 and proceeding exactly 
as did Goldstein and Siegel the hodograph and complex potential planes 
are as shown in Figure 3-1 end 3-2. The function which maps the 
rectangular region T in the T-plane (figure 3-3) into the hodograph 
plane is 

C (T) = V q e l(TT/2 “ a) fi(T;§ H ,k) Vfi(T;? G> k) [0(T;§ A> k)] e/TT (1) 

where the Q function is as defined in Appendix A. Taking the argument 
of this function and requiring that the direction of the flow along the 
plate must equal the argument, it can be shown that 
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Physical flow field, 2-d tilted jets, 
z = x + iy 
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The remainder of the analysis proceeds exactly as before and the 
final solution for displacement in the flow field becomes 



2 = /T\ ia k'_ |- sn (ri F’ kl) + 6 R /6 Ll 

A VV J e n L 1 + 6„/6, j 

o R L 



X 



T 



/ 

V iK ’ 



(3) 

(k sn snT-1) (1-k snT) dT 
*vh 

dnT[l-dn(n F ,k') snT] n(T;^,k) M T ;§ G ,k) [fi(T ;§ A ,k)] P/lT 



There are in this case six parameters which determine the flow field: 
plate height, plate width, plate angle, nozzle height, nozzle spacing and 
nozzle angle. 

The plate and nozzle angles appear in the solution but the other 
parameters must be calculated. Following the example of Goldstein and 
Siegel the solution is determined by integrating along the sides of the 
rectangular region T in the counterclockwise direction. Thus the distance 
between two points in the physical plane is the difference between the 
line integral from the point G in the T-plane (origin in the physical 
plane) to each of the points in question. If all distances are non- 
dimens iona lized by the nozzle width the solution can be obtained using 
the known plate and nozzle angles and the four following input parameters: 
k, 6 T /6 n , ru/K 1 . and £./K. These parameters all have values ranging 
between zero and one, but no representative values were given in Ref. 3. 
Upon direct communication with the Lewis Research Center it was learned 
that no values of these parameters and their corresponding physical con- 
figurations were available. 
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IV. COMPUTER ANALYSIS 



Using the working formulas in the appendixes of Refs. 3 and 4, two 
computer programs were written. The first program calculated the plate 
height and width and nozzle height and spacing for a given set of input 
parameters. The second program was merely a variation and extension of 
the first. For a given set of input parameters and experimental data the 
second program plotted the theoretical flow field, the pressure coeffici- 
ents on the ground and plate and the velocity profile in the nozzles along 
with the data points. A copy of the second program follows Appendix A. 

The plate and nozzle angles were read from data cards. In the second 
program the input parameters were also read from cards, but in the first 
program the parameters were initialized and iterated in steps of one 
tenth. The complete elliptic integral of the first kind, the complimentary 
modulus and the complimentary integral were all calculated from k. The 
various quantities and functions were then evaluated as indicated in the 
program. 

The modified equation for § is equation (2) of this thesis. The 
modified equation for 0 is 



® (i) (ti) = \ - a + Q^Lti) + \ Qq 1) (ti) + ^ Q^ 1} (ti) i = 3,4 



(4) 



These modifications provide for solutions with the nozzles tilted. The 



>(i), 



(3) 



/f s are angles defined in Ref. 4. It should be pointed out that 
must be taken as a negative angle in the third or fourth quadrant, while 

(4) 

must be taken as a positive angle in the first or second quadrant. 
The height and width of the plate were calculated directly by inte- 
gration. The nozzle spacing was found by first locating the poir, B at 
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the end of the plate and then adding to that the horizontal separation of 
the points B and C. The nozzle height was found similarly by adding the 
vertical separation of the points D and E to the width of the flow to the 
right . 

At this point the values calculated were printed out. In the first 
program the parameters were iterated here and a new configuration was 
calculated. In the second program plotting was begun at this point. 

The data points were plotted first and then the theoretical free 
streamlines were plotted. A similar procedure was followed for the 
pressure coefficients and the nozzle exit velocity. The velocity was 
made negative to fix the direction downward, so that the plotted profile 
was as if it were being viewed from in front of the aircraft. 

All integrations were performed by the Gauss nine-uoint quadrature 
formula. When infinite series were specified only the first five terms 
were taken. In the working formulas presented in Ref. 3 quite often 
square roots and odd powers were taken of the negatives of numbers. On 
the intervals indicated this would be quite proper since the numbers 
would be negative values before having the sign changed. Since the 
programs were to be used to find possible combinations of the input 
parameters, quite often, if the parameters were wrong, the result would 
be an attempt to take the logarithm or square root of a negative number. 
To prevent the computer from ending the calculation with error messages, 
absolute values were taken of all functions under radicals. This was 
equivalent to the operation intended in the original notation, but 
avoided the problems in calculation. 

In the first program statements were included in all the functions 
to notify the operator if a denominator was zero. The second program 
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has retained the check on the denominator and sets the function to a 
nominal value of one thousand in the event the denominators are very 
small. This is intended only to avoid underflow errors at the limits of 
the integration; no notification is given to the operator in this case. 

The second program is therefore intended to be used only with combinations 
of the parameters which are known to satisfy all the equations. 
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V. DESIGN CONSIDERATIONS 



A model of a V/STOL aircraft using two-dimensional slot jets in an 
inviscid, incompressible fluid was required to verify the theoretical 
results. The test apparatus was designed to approximate these conditions 
as closely as possible. 

Air at low pressure was used to represent the incompressible, inviscid 
fluid. The two-dimensional slot jets were produced by two rectangular 
nozzles with a cross-sectional width to length ratio of four to one. By 
confining all measurements to a plane through the centers of the two 
nozzles, end effects could be ignored. 

To avoid entrainment problems the theory required that the nozzle 
half-spacing and height be only a few jet widths. This condition 
suggested the use of a large nozzle width so that there would be enough 
room to probe the flow field beneath and between the nozzles. Also to 
permit a significant number of data points for plates with half widths 
as small as one fourth of a jet width, a large nozzle width was required. 
The test apparatus was designed with nozzles two inches wide. 

Even with low nozzle exit velocities a large flow rate was necessary 
to supply the system. Thus, the nozzles were connected to the laboratory 
turbine compressor. Each nozzle was attached to the air supply through a 
flow straightening, diffuser chamber. Rotameters preceded these chambers 
in the air lines. Thus, air flow conditions could be monitored easily 
throughout the test period. As air entered the diffuser chamber it was 
deflected by three layers of wire screen. The air then passed through a 
layer of honeycomb material into the nozzle, thus producing a uniform 
velocity distribution throughout the nozzle. 
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Laboratory model of V/STOL aircraft using 
symmetrical two-dimensional slot jets 

FIGURE 4 




Flow straightening chamber 
FIGURE 5 
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The flow straightening sections also served as a mounting for the 
.nozzles. If desired, the nozzle shape and size could be changed quite 
easily. 

A mounting arrangement was designed to not only allow the nozzles to 
be moved vertically and horizontally but also to be rotated in the vertical 
plane, thereby representing a tilted jet configuration. This freedom of 
movement was accomplished by connecting the flow straightening section to 
a vertical rod through a swivel which could be locked at any angle. The 
vertical rod passed freely through a metal block, and it too could be 
locked in any position. The block moved horizontally along parallel bars 
and could be held firmly in position with a set screw. The parallel bars 
were mounted in a rigid framework which was supported by four legs. 

Directly underneath the framework and the nozzles was a large plexi- 
glass plate which represented the ground. Smaller plates representing 
fuselages of various sizes and shapes were made to be interchangeable on 
a bracket which was suspended from the framework between the nozzles. 

The vertical position of the plate bracket was varied with a threaded 
traverse. The ground plate and the fuselage plates contained a row of 
pressure taps across the center section. For future studies with smaller 
nozzles and higher velocities, significant pressure measurements would be 
made possible, but for these low velocity tests the primary measurements 
were made with a hot-wire anemometer. 

To support the hot-wire probe a traversing mechanism was mounted on 
the apparatus. The mechanism allowed the probe to be held at a constant 
depth into the flow while traveling both horizontally and vertically 
throughout the field. 
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Nozzle mounting arrangement 
FIGURE 7 




Traversing mechanism 
FIGURE 8 
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VI. PROCEDURE 



A list of plate heights and widths and nozzle heights and spacings 
for which the corresponding input parameters were known was compiled. 

From this list a few representative configurations were chosen for 
experimental verification. 

The plate size, height and angle and nozzle spacing, height and angle 
were all set and recorded. The compressor was then started and an equal 
air flow was routed to each nozzle through the flow meters. A straight 
hot-wire probe was placed under one edge of one of the nozzles. The 
voltage and location was recorded. Then by following the motion of a 
fine thread placed in the flow field, the probe was moved horizontally 
in convenient intervals and then moved vertically until the same voltage 
was again observed. Since only velocity magnitudes were of interest, the 
probe was always rotated to obtain a maximum reading. By following this 
procedure from both sides of a nozzle and from the edge of the upper plate 
the coordinates of the lines of constant voltage (free streamlines) could 
be recorded. 

The probe was then moved as close to the upper plate as possible and 
oriented to read the velocity in the direction of the plate. By moving 
the probe- in convenient intervals along the plate the velocity distri- 
bution along the plate could be recorded. Similarily the probe was moved 
along the lower plate and the voltage and horizontal position were 

>v : 

recorded . 

This procedure was followed with each configuration. Tests were run 
with Reynolds numbers varying from 1600 to 7000. Simultaneously pressure 
measurements were made using a micromanometer capable of indicating 
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pressures as low as .001 inches of water. The pressure and location were 
recorded for comparison with the hot-wire data. 

The data collected in this manner was sufficient to plot the free 
streamlines of the flow field, the nozzle exit velocity profile and the 
pressure coefficient along the plate and ground. The input parameters for 
each configuration were read into the program along with the experimental 
data. The experimental results for each configuration were then plotted 
along with the corresponding theoretical results. 
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VII. RESULTS 



The results of the experiments are presented in the following pages. 
Table I is a list of some of the input parameters for which the corres- 
ponding configurations were found. 

Figures 8-1 through 12-1 are typical free streamline plots for various 
configurations. Hot-wire data points are represented by triangles. The 
coordinates are nondimensionalized with the nozzle width. 

The pressure co fficient along the ground is shown in Figures 8-2 
to 12-2. Only the abscissa is nondimensionalized. The pluses represent 
data points calculated from manometer readings. 

The pressure coefficient on the plate is shown in Figures 8-3 through 
12-3. Again, only the abscissc is nondimensionalized. 

Figures 8-4 through 12-4 are the plots of the nozzle exit velocity 
profile. The velocity is nondimensionalized with respect to the maximum 
velocity recorded for that test. The abscissa is one nozzle width in 
length. The velocity is plotted negatively so that the plot represents 
the actual direction of the flow. 

Figures 13 and 14 clearly indicate the nature of the flow field 
beneath the nozzles. 
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TABLE I 

Input parameters and corresponding configurations 
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FIGURE 8-1 Free streamline plot 
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FIGURE 8-2 Pressure coefficient along the ground 
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FIGURE 8-3 Pressure coefficient along the fuselage 
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FIGURE 8-4 Velocity profile in the nozzle exit 
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FIGURE 9-1 Free streamline plot 
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FIGURE 9-2 Pressure coefficient along the ground 
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FIGURE 9-3 Pressure coefficient along the fuselage 
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FIGURE 9-4 Velocity profile in tlie nozzle exit 
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FIGURE 10-1 Free streamline plot 
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FIGURE 10-2 Pressure coefficient along the ground 
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FIGURE 10-3 Pressure coefficient along the fuselage 
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FIGURE 10-4 Velocity profile in the nozzle exit 
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FIGURE 11-1 Free streamline plot 
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FIGURE 11-2 Pressure coefficient along the ground 
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FIGURE 11-3 Pressure coefficient along the fuselage 
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FIGURE 11-4 Velocity profile in the nozzle exit 
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FIGURE 12-1 Free streamline plot 
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FIGURE 12-2 Pressure coefficient along the ground 
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FIGURE 12-3 Pressure coefficient along the fuselage 
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FIGURE 12-4 Velocity profile in the nozzle exit 
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Tufts of yarn across nozzle exit 
FIGURE 13 




Strips of paper attached to ground plate 
and nozzle 

FIGURE 14 
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VIII. DISCUSSION OF RESULTS 



Many combinations of the input parameters produced impossible con- 
figurations (negative plate heights, etc.)* In fact there were far more 
extraneous solutions than useful ones. There were no combinations found 
which produced a configuration wi h the nozzle exit lower than the plate. 
Very few reasonable configurations resulted from the maximum or minimum 
values for the parameters. 

Most of the configurations of interest were found to have input 
parameters in the following ranges : 

.3 < k < .6 

•'< 6 l /6 R^ ' 5 
.6 < rip/K' < .9 
.1 < § a /K < .9 

There were possible configurations outside of these ranges, but most of 

them had extremely small plate widths. Variation of any one input 

* 

parameter was not found to change any single physical parameter. 

Decreasing the plate angle had a tendency to decrease the plate width 
and increase the plate height. Tilting the nozzles caused an increase in 
plate and nozzle height but decreased plate width and nozzle spacing for 
a given set of abstract input parameters. Neither of these variations had 
any outstanding effect on the theoretical pressure coefficients. 

The experimental data for the streamlines agree quite well with the 
theoretical plot. The data points between the plate and nozzle often 
seem to follow a different streamline than the theory. This is a result 
of the fact that the test plates were of finite thickness, while the 
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theoretical plates were infinitely thin. Figure 9-1 shows entrainment 
effects for very high nozzle positions. In some cases the theoretical 
plot can be very impractical. Notice in Figure 11-1 that the streamline 
from the edge of the plate turns and passes through the plate. 

In most cases the theoretical pressure coefficients along the ground 
and plate were quite close to the actual situation. The micromanometer 
data points are somewhat questionable. The pressures were extremely low 
even for the instrument which was used. The theoretical coefficient is 
assumed to approach zero after the last point plotted, as the measurements 
indicate . 

The nozzle exit velocity profiles for very low velocities were 
approximately the same shape as those from the theory, but as the exit 
velocities become higher the profile becomes much flatter. In all cases 
it should be remembered that the velocity data points are normalized with 
the maximum velocity for that particular test. In most cases this maximum 
velocity was found to be near the ground two or three nozzle widths from 
the high pressure region beneath the fuselage and nozzle exit. Hence, the 
data points tend to be displaced from the theoretical profile which assumes 
that the maximum velocity is reached at the edges of the nozzle exit. In 
all cases the profile was somewhat affected by the presence of the wall of 
the nozzle. 

The motion of the tufts of yarn placed in the nozzle exit (Figure 13) 
show the manner in which the flow divides to the left and right of the 
nozzle. The strips of paper attached to the nozzle exit and ground plate 
in Figure 14 show the stream of flow away from the nozzle an the upward 
air currents under the fuselage. 
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IX. CONCLUSIONS 



The experiments have not only verified the theory, but have also 
pointed out its limitations. The analytical results were quite accurate 
for nozzles very close to the plate and the ground. As the nozzles were 
moved higher the effects of entrainment became more noticeable. The 
pressure distribution along the ground was not as uniform as that pre- 
dicted by the analysis. 

The greatest limitation on the theory, however, is not the restricted 
nozzle positions, but the fact that the analytical solution depends upon 
four parameters in an abstract plane. Since there is no direct connection 
between any one of these parameters and the physical parameters of plate 
height and width and nozzle height and spacing, the theoretical results 
are only possible after an iteration process to find the proper input 
parameters. A vast amount of computer time is required to solve for 
relatively few configurations. 
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X. RECOMMENDATIONS 



The theoretical solution could possibly be extended to include other 
fuselage shapes. However , perhaps a more practical approach would be to 
develop analytical solutions in terms of physical parameters. A solution 
of this type would be a useful tool to the aircraft designer. 

Further tests should be carried out with higher nozzle exit velocities 
to determine the effects of extremely turbulent flow underneath the air- 
craft. Also of interest would be the effect of crosswinds on the flow and 
the pressure distribution. 

Reference 1 mentions the use of a special takeoff grating to virtually 
eliminate entrainment effects on aircraft using a single central jet 
exhaust. No information was available concerning the effects of such a 
grating on multi-jet configurations. Tests of a takeoff grating mounted 
on the existing apparatus would be. of interest to determine the resulting 
changes in ground effect. 

There were no reasons given for restricting 6 T /6 D to values less than 

L R 

one. Obviously there must be solutions for nozzle exits lower than the 
plate height. If the equations can be satisfied by values of 6 ^/ 6 ^ 
greater than one, perhaps this will lead to the lower nozzle positions. 
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APPENDIX A 



To obtain analytical solutions for the flow field beneath a V/STOL 
aircraft Siegel and Goldstein [Ref. 3] considered a simplified two- 
dimensional model of the flow. The following is a condensation of their 
analysis . 

The aircraft under consideration has a fan pod or jet exhaust in each 
wing. As the two jets strike the ground they divide into flows to the 
left and right. The colliding flows beneath the aircraft rise and 
strike the fuselage. The bottom of the fuselage is represented by a 
flat plate which can be at an angle with the ground. Because of symmetry 
only half of the flow field is considered. 

The flow is assumed to be incompressible, inviscid, irrotational , 
isothermal and steady. The velocity at the nozzle exit is taken to be 
perpendicular to the exit plane. The surrounding fluid is at rest and 
entrainment is ignored. If the nozzle exits are within a few jet widths 
of each other and of the ground, the inviscid solution with no entrain- 
ment is a reasonable approximation of the actual flow conditions. With 
these restrictions the flow has been reduced to a free streamline prob- 
lem with constant pressure surroundings. 

Figure Al-1 shows the physical boundaries of the flow field. Defining 

the complex conjugate velocity as 

Q = u - iv 

and the complex potential as 

W = 5 + ity 

the flow regions in the hodograph and complex potential planes are shown 
in Figures Al-2 and Al-3. 
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FIGURE A 1-2 Hodograph plane 



FIGURE Al-3 Complex potential plane 
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The solution in the physical plane is obtained by performing the 



integration 



z 




dW + constant 



(A3) 



In order to perform this integration Q and W must be expressed in 
terms of the same variables. This is accomplished by mapping the 
rectangular region T of the T-plane onto the Q and W planes. 

The function which maps the region F onto the hodograph plane is 

C (T) = V o e 3 CV2) fi(T;§ H ,k) Vo(T;§ G> k) [Q(T ;§ A> k) ] ?/tT (A4) 

The Omega functions are defined in Ref. 4 to be 



ft(T ;§ r >k) = 



^(fit <T+e r -2K)l|f9 



(A5) 



where 9^ is the fourth theta function which is defined in Ref. 8 as 
follows : 



0, (z | T) = 1+2 S (-l) n 



im 






cos 2nz 



n=l 



(A6) 



The relationship between § , § and § must satisfy the requirement 

AG H 

that on the plate the argument of £(T) must equal the direction along the 
plate. By using the formulas for the arguments of the Omega functions 
presented in Ref. 4 it can be shown that 



V4 + n 5 H ^ n h 



(A7) 



By applying a Schw.rz -Chris toffe 1 transformation the upper half of 
an intermediate t-plane (not illustrated) is mapped onto the complex 
potential plane. A second transformation maps the T-plane onto the 
upper half of the t-plane. By eliminating t between these two trans- 
formations, the function which maps the T-plane onto the W-plane is 
found to be 



54 



dW 

dT 



i.A 



(A8) 



(k sn § sn T - 1) dn T 

(1 + k sn T) [ l-dn(r]„,k 1 ) sn T] 

r 

where A is a constant . 

Using partial fractions, integrating equation (A8) and applying the 
condition that the flow through the nozzle exit is equal to the differ- 
ence in the stream function across the exit, the value of A can be found 
in terms of the average nozzle exit velocity. 

Noting that the flow to the left beneath the nozzle is equal to the 
change in the imaginary part of the complex potential at F and requiring 
the velocity at that point to be the free streamline velocity we have 



_ k + dn Cr| F ,k' ) - k (1 + sn § R ) 

6 L V o I A k + dn (ri F , k ' ) - k (1 + sn [ 1 - sn (r| F , k ' ) ] 

Since continuity requires that 

6 d V + 6 T V = AV T (A 10) 

R o L o I 

we can eliminate 6 T from equation (A9) and have 

Xj 



6 r v i 



k (1 + sn § H ) sn (r| F ,k') 



A V q k + dn (ti f , k ‘ ) - k (1 + sn § R ) [ 1 - sn (r) F ,k ' ) ] 
By dividing equation (All) by equation (A9) and rearranging 

6 r /6 l dn <dV k '> ' k sn (’Hp ’k ' ) 



(All) 



Sn ? H = 



6 r /&L k + k sn (ri F > k ’ ) 



(A 12) 



Since Q and W are both known in the T-plane the solution of equation 
(A3) becomes 



1 dW 



z = 



C(T) dT 



dT 



(A 13) 



§ G +i Rl 
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or after substitution 





psn (r| F ,k') + 6 r /6 l -j 


VV J TT 


L 1 + 6 / 6 t J 



(A 14) 



(k sn § H sn T - 1) (1 - k sn T) dT 



§ G + iKl 



dn T [ 1-dn (r| p ,k ' ) snT] 0(T;§ H ,k)7o(T;§ G ,k) [Q(T ;§ A ,k) ] 0/rr 



Using this relation the ratio of the average velocity in the nozzle 
exit to the free stream velocity is 



- = M [• 



1 + 6 R /6 L 



sn (r| F , k ' ) + 6 r / 6 l - 



-] 



(A15) 



where 



M 



K 



(1-k sn ? H sn§) (1-k sn §) d? 



]" J dn£ [ l-dn(n F ,k- )sn? ]Q(§ ;§ H ,k)Vfi(§ ;§ G> k)'[fi(§ ;? A ,k) ] 0/rr 

" K 



(A 16) 



The 



pressure coefficient at a point (x^, y ) is 



C = 



p( v y{> p o 

1/2 p V 2 



(A 17) 



Using Bernoulli's equation the pressure coefficients on the ground 
and on the fuselage are 

• , \~|2 



C = 1 
P>8 






-K < § < §, 



and 



(A 18) 



5 A <?< 



K 



The flow field is governed by five physical parameters: plate height, 

plate width, plate angle, nozzle height and nozzle spacing. Only the 
plate angle 3 appears in the solution. Hence the other parameters must 
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be found in terms of four parameters in the T-plane. Siegel and Goldstein 
used k, r|p/K T and E^/K as input parameters for a computer program 

to plot nume Leal results. All four parameters have values between zero 
and one but no representative values and their corresponding physical 
configurations were available. 
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EQN 3 5288 (MODIFIED FOR NOZZLE TILTING) 

XIG=EL + ( TEL/ PI )*< ALPHA + 3ETA/2) ) - 2o *XI H- 2 , -5= BETA* X I A/P I 

BIGM IS INTEGRAND OF EQN 18 5288 

CALL QT 9B( -EL , EL , B I GM , R4 , E L , ELP , EMOD , SNH , X IH , X I A , X I G , DNF , BET A , P I 
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R = R+ol561735*(FCT( A+C , AC H, T , EL , SMQ, EMOD , SNH, X IH,X I A , XIG, ETAF , BETA 
1, ALPHA )+FCT ( A-C, ACH , T , EL , SMQ , EMOD , SNH , X I H , XI A , X I G , ET AF , B ET A , ALPH A 
2 ) ) 

R=6* (R+o 16511S7*FCT(A, ACH,T,EL ,SMQ,EMOD,SNH,XIH,XI A,XIG, ETAF, BETA 
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